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Experimental data on sulphur and oxygen nuclei interactions with pho- 
O ■ toemulsion nuclei at the energies of 200 and 60 GeV/nucleon are analyzed 

with the help of a continuous wavelet transform. Irregularities in pseudo- 
rapidity distributions of narrow groups of the secondary shower particles 
in the mentioned interactions are observed at application of the second 
order derivative of Gaussian as a wavelet. The irregularities can be in- 
terpreted as an existence of the preference emission angles of groups of 
particles. Such an effect is expected at emission of Cherenkov gluons in 
nucleus-nucleus collisions. Some of the positions of the observed peculiar- 
^ ! ities on the pseudorapidity axis coincide with those found by I.M.Dremin 

et al. (I.M.Dremin et al. Phys. Lett., 2001, v. B499, p. 97). ) 

In the last decade a new type of mathematical analysis of data, the so-called wavelet 
analysis has become very popular in various branches of science and engineering [0-0. 
Mainly, it is applied for analysis of time series coming from geophysics, meteorology, 
. astrophysics, and so on (applications in aviation, medicine and biology see in ||). 

Qh! The wavelet decomposition or wavelet transform of a function f(x) is its decomposition 

J|h I on an orthogonal functional family of special form ||; 

Wv(a, b)f = — ^= r f(x)V a , b (x)dx, (1) 
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where \I/ is called a wavelet, b - a translation parameter, a - a dilation parameter or a 
scale, and C$ - a normalizing constant 
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where 4f(u) - Fourier transform of ty(x). 

The derivatives of the Gaussian function are often used as wavelets 

Jn 

v[/(x)^^(x) = (-ir +1 ^e-^ 2 , n>0, C 9n = 2n(n - 
The first two wavelets are well known: 
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The second one is called Mexican Hat wavelet (MHAT). 

As seen, the wavelet transform puts in correspondence to the function of one variable, 
f(x), the function of two variables, W\y(a, b). Until recently, a presentation and an analysis 
of a function of two variables was quite a difficult job, and only modern computers with 
their 3D-graphics allowed one to implement completely the method of wavelet analysis. 

There is a discrete analogy of the continuous wavelet transform (see ||). It was used 
in elementary particle physics for the study of events of cosmic rays interactions with 
materials - at the analysis of particle pseudorapidity distributions 0. The wavelet coef- 
ficients (Wy(a, b)) of energy distributions predicted by different models of multi-particle 
production were studied in || . Possibilities of the wavelet analysis in searching for mani- 
festation of the disoriented chiral condensate in the pseudorapidity distributions of neutral 
particles fraction were considered in M. 



In papers JT0[ the wavelet transform was used for pattern recognition in Pb + Pb- 
interactions at the energy of 158 GeV/nucleon. Structures were found in the angular 
distributions of secondary particles that can be interpreted as irradiation of Cherenkov 
gluons. The last publications in this research are devoted to experimental search for the 
disoriented chiral condensate in nucleus- nucleus interactions [ |TT|j . 

Most of the mentioned papers suffer from the apparent lack of quantitative results 
which caused either by uniqueness of nature phenomenon or low statistics of analyzed 
data. It is connected partly with specific properties of the wavelet analysis itself (see 



12 1) . A regular method of wavelet analysis application in particle physics is needed. Our 
paper presents experience of using the wavelet transform for analysis of more than 2000 
interactions of nuclei with nuclei at high energies. 

Experimental data were obtained at horizontal irradiation of NIKFI BR-2 nuclear 
photoemulsion by sulfur and oxygen nucleus with the energies of 200 and 60 GeV per 
nucleon at the CERN SPS. The sensitivity of emulsion was about 30 grains per unit 
length of 100 \xm for single charged particles with minimal ionization. 

Primary interactions were found by along-the track double scanning: fast in the for- 
ward direction and slow in the backward direction. Fast scanning was made with a 
velocity excluding any discrimination of events in the number of heavily ionizing tracks, 
slow scanning was carried out to find events, if any, with little changed and unbiased pro- 
jectile nucleus. Upon rejecting events of electromagnetic dissociation and purely elastic 
scattering in the total sample, 884 events of 5* + Em and 504 events of O + Em at the en- 
ergy of 200 GeV/nucleon, and 884 O + Em interactions at the energy of 60 GeV/nucleon 
were selected for a further analysis. In each event the polar angles 9 and azimuthal angles 
ip were measured. 

Shower particles, or the so-called s-particles - single charged particles with a velocity 
P = v/c> 0.7, are considered in this study. The s-particles consist mainly of produced 
particles and single charged nuclear fragments. A special separation of single charged 
fragments was not done. 

According to Ref . |I(J , a distribution of the secondary particles on the pseudorapidity 



T] = — ln(tan(0/2)) in an event was presented as 
where N is multiplicity of s-particles in the event, and rji is pseudorapidity of i-th particle. 
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A wavelet transform of the function (fj) gives 

So, the function of two variables is brought into correspondence with each particle. 
Wavelet spectra of the event (W$(a, b)) is a sum of such functions. 

As an example, Fig. 1 presents wavelet spectra of event with 6 particles having 
pseudorapidities r\\ = 1, r/ 2 = 2.75, r/ 3 = 3.25, ^4 = 5, r/ 5 = 6, r/ 6 = 7. The values of 
translation parameter b are put on the X-axis, dilation parameter, a - on Y-axis. The 
values of the wavelet coefficients are depicted in a gray-level scale: the high values of 
the coefficients are of light shade while the lower ones are darker. As seen, in #2 wavelet 
spectra at the scale lower than 0.3 all particles are distinguished. At the scale larger than 
0.5, the particles 2 and 3 can not be resolved. At a > 1 particles 4, 5, 6 can not be 
resolved. At a > 2 one could expect a fusion of particles 1,2, and 3. However, this does 
not take place due to small yield of particle 1 in the wavelet spectra at a large scale. So, 
the wavelet transform of g% allows one to study the particle clusterization. 
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Figure 1: Wavelet spectra, energy densities, and scalogramms 
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Let us note that the positions of the local maximums of Wy(a, b) (a ~ 0.5 and 6 = 3, 
a ~ 1 and b = 6) on the b-axis are connected with positions of the centers of the groups 
of particles 2, 3 and 4, 5, 6. The corresponding scales reflect the width of the groups 
on the pseudorapidity axis. Positions of the local minimums give the centers of the 
pseudorapidity splits between the groups. 

There are a lot of wavelets. Different representations of the 6 particles of the test 
example are shown in Fig. 1. There are also energy densities (Wq,(a, b) 2 ) and scalogramms 
E w (a) = J W^,(a,b) 2 db, which are often used in practice. 

The gi wavelet (see (2)) has a minimum at a negative value of its argument and 
maximum - at positive value. Positions of these extreme points at different scales can be 
seen quite well in the energy density spectra. The g 2 wavelet has 3 extreme points. All of 
them are presented in the energy spectra (see the region of the 1 particle). Scalogramms 
give generalized imagination of the extreme points. 

As seen in the figure, different wavelets give different presentations of particles. The 
more "usual" ones are obtained at the application of even wavelets. Perhaps, odd wavelets 
can be useful at automatic data processing. 

At the first glance, the g\ wavelet allows one to locate the regions of unhomogeneouslity 
in the particle distribution. The positions of the local maxima of the energy spectra on 
the b-axis (6 = 2.5 and b = 3.5) mark the group of the particles 2 and 3, and the positions 
on the a-axis show the group size. Though, the following local maxima are disposed 
at a = 4 and b ~ 0, b ~ 9. There is no selection of the group of the particles 4, 5, 
6. At the same time, there are irregularities (extreme points and inflection points) in 
the corresponding scalogramm which are probably connected with characteristics of the 
particle group. 7 There are extreme points in the scalogramm obtained with help of the 
g 2 wavelet at b ~ 0.5 and 1.2 associated with the group characteristics. So, we believe 
that scalogramms can be used for fast search of the particle group. 

The g 4 wavelet has analogous properties. Though, the characteristics of the corre- 
sponding scalogramms are not so strongly correlated with the particle characteristics. 
Thus, below we will use mainly the g 2 wavelet. 

We started our analysis with study of the particle distribution on the pseudorapidities 
in all interactions. The histogrammed distribution for S + Em interactions at the energy 
of 200 GeV/nucleon is presented in Fig. 2. There are also g 2 and g± wavelet spectra 
at different scales. It should be noted that the wavelet transform was applied to raw 
experimental data, not to the histogrammed distribution. Thus, a fine structure of the 
spectra can be observed at small scales. The structure becomes more regular at large 
scales, and turns to the wavelet function in the limit of large a. At a ~ 0.4 one can see 
3 maxima in the g 2 spectra at rj ~ 2, 4 and 8. It seems that the maximum at rj ~ 2 is 
connected with the particle production in the target fragmentation region; the maximum 
at rj ~ 4 - with the particle production in the central region, and the maximum at rj ~ 8 
- with spectator fragments of the projectile nucleus. In reality, the three maxima only 
reflect what can be seen with the naked eye - the left wing of the distribution is different 
from the right one. One can clearly see a change of the slope on the right wing at rj ~ 6. 
At smaller scales the structure is more rich and does not allow such simple interpretation. 
In order to find selected scales in the interactions and to analyze a fine structure of the 
events, we turned to a study of the energy density and scalogramm. 

The spectrum of the energy density reflects peculiarities of W^(a, b) function and can 
be used for searching the particle group. Though, definition of maxima and minima of the 
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function is a quite complicated task. So, at the first stage we concentrated our attention 
on the energy distribution on scales, on scalogramm because it is a ID-function. 





Figure 2: Wavelet analysis of s-particles pseudorapidity distribution in S + Em interac- 
tions at the energy 200 GeV/nucleon 



Scalogramms reflect characteristic features of events according to the test example. 
For example, the scalogramm of the g-i spectrum has the minimum at a ~ 1.1 associated 
with an average distance between the groups of the particles 2, 3 and 4, 5, 6, and the 
maximum at a ~ 0.5, connected with the most compact group of particles 2, 3. It is 
possible to write out an analytical expression for the scalogramm using the g 2 wavelet 
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and the distribution of the form (fj). 
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It is seen that the scalogramms is a statistical averaged squared distance between particles. 
Thus, we expected that an analysis of the scalogramms allowed us to find characteristic 
scales. 

Distributions of local minimum and maximum positions in the scalogramms of the 
events of S + Em interactions at the energy 200 GeV/nucleon are given in Fig. 3. Though 
there are some irregularities in the distributions, we can not say that they have a statisti- 
cally guaranteed meaning. In addition, the analysis of connection between a m i n and a max , 
and characteristics of the real events did not show any regularity. However, we believe 
that the study of a larger volume of experimental data will produce interesting results. 
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Figure 3: Distributions of the extreme points of the scalogramms of the S + Em events 



The next step of our study was search for extremum points of function Wy(a, b) of the 
real events. Distributions of points on the dilation scale in our interactions are presented 
in Fig. 4. At the first glance, the distributions have no any characteristic peculiarities 
- bumps or pits. We only can mark that the distributions can not be described by a 
simple exponential function. It is necessary to use at least 2 exponents at fitting the 



distributions. Simulation of pp- interactions with the help of the HIJING program [13 



has shown that the appearance of the exponent with a larger slope is connected with 



6 



production of jets of particles. So, we consider our distributions as a manifestation of jet 
existence in nucleus-nucleus interactions. 




0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 4: a max distributions in the S + Em and + Em interactions at the energy of 
200 GeV/nucleon, and in the O + Em interactions at the energy of 60 GeV/nucleon 
(histograms 1, 2 and 3, respectively). The distributions 1, 2, and 3 are multiplied by 10 2 , 
10 1 and 10°, respectively. 



Research in the distributions of local maximum of Wy(a, b) on b gives more interesting 
results. The Fig. 4 shows the distributions of all our interactions at a max > 0.05 (1), 
a-max > 0.1 (2), a max > 0.2 (3), and a max > 0.3 (4). As seen, the peculiarities of the 
distributions at r\ ~ 1.5, r\ ~ 2, r\ ~ 3, r\ ~ 3.5, r\ ~ 5 are located at the same positions 
for different interactions. The peculiarities, as it can be seen, are connected with narrow 
groups, where a max < 0.05. We interpret them as an existence of preference emission 
angles of the groups of the particles. Let us note that the positions of some irregularities 



found by us coincide with those observed early in Ref. [[10] at the study of Pb + Pb 
interactions at the energy of 158 GeV/nucleon. Unlike Ref. [|Hj] where 5 mostly central 
events were analyzed, we present the results for more than 2000 events. 
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Research on narrow groups of the particles with help of traditional methods was the 
third step. Fig. 5 shows a distribution on a pseudorapidity interval between neighboring 
particles in an event. It was observed that there were pairs of the particles with close 
pseudorapidities, A < 10~ 5 ! Their number increases statistical fluctuations. The distri- 
bution of the centers of the pairs on r\ has two bumps at r\ ~ 3 and 4. So, we can conclude 
that the peculiarities of the wavelet spectra observed by us at small scales are connected 
with the narrow groups of the particles. 

One can suppose that such pairs occur due to data input error when two or three 
entries in the event record are corresponded to the same particle. In this case, such ghost 
"particles" must have identical r] and tp. The fraction of these pairs is lower than 20 % 
among all the observed narrow pairs. The other pairs have different values of azimuthal 
angles. The distribution of the particles of the narrow groups on tp (see Fig. 6) has no 
clear peculiarities and reflects the methodical drawback of the photoemulsion experiments 
- a poor identification of particles flying perpendicularly to the emulsion plate (at tp ~ 0° 
and 180°). The distribution of the azimuthal angle difference of the particles of the narrow 
group has no irregularities either. So, these pairs can not belong to a jet of particles. The 
nature of such pairs is not clear for us. 

Summary 

1. The continuous wavelet transform has been used for the analysis of more than 
2000 events of nucleus- nucleus interactions at high energies (S + Em and O + Em 
interactions at the energies of 200 and 60 GeV/nucleon). 

2. It is shown that the maxima of W^,(a,b) obtained with the help of the gi wavelet 
are associated with the groups of particles. 

3. It has been found that the distribution of the group of particles on scales in the 
interactions under the study is geterogeneous what can be caused by jet production. 

4. It is observed that the distributions of the groups on pseudorapidities have irregu- 
larities, there are preferences of emission angles of the groups. 

5. The pairs of particles with close pseudorapidities, A < 10~ 5 , are found for the first 
time. 

The nature of the peculiarities observed by us is not clear yet. 

The authors are grateful to the members of the EMU-01 collaboration for their kind 
permission to employ the experimental data analyzed in this paper. One of the authors 
(V.V.U.) thanks RFBR (grand N 00-01-00307) and INTAS (grand N 00-00366) for their 
financial support. 
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Figure 5: b max distributions in the S + Em and + Em interactions at the energy of 200 
GeV/nucleon, and in the + Em interactions at the energy of 60 GeV/nucleon. Ab is 
step of the histogramming. The distributions 1-4 are obtained at a max > 0, a max > 0.05, 
a ma x > 0.1 and a max > 0.2, respectively. 
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Figure 6: Distribution on the pseudorapidity interval between the neighbouring particles 
in all studied interactions. 
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Figure 7: a) The azimuthal angle distribution of particles with close pseudorapidities. 
The distribution of all particles is presented by solid line, b) The distribution on the 
azimuthal angles difference of particles with close pseudorapidities. 
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